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Abstract 

Gamma-rays are, with antiprotons, a very efficient way to derive upper limits on 
the density of evaporating black holes. They have been successfully used in the 
last decades to severely constrain the amount of Primordial Black Holes (PBHs) 
in our Universe. This article suggests a little refinement, based on the expected 
background, to improve this limit by a factor of three. The resulting value is : 
VLpbh < 3.3 x KT 9 . 
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Small black holes could have formed in the early Universe if the density con- 
trast was high enough (typically 5 > 0.3 — 0.7 , depending on models). Since 
it was discovered by Hawking [1] that they should evaporate with a black-body 
like spectrum of temperature T = he 3 /(8irkGM), the emitted cosmic rays have 
been considered as the natural way, if any, to detect them. Those with initial 
masses smaller than »5x 10 14 g should have finished their evaporation by 
now whereas those with masses greater than a few times M* do emit nothing 
but extremely low energy massless fields. The emission spectrum for particles of 
energy Q per unit of time t is, for each degree of freedom, given by : 



where k is the surface gravity, s is the spin of the emitted species and T s is the 
absorption probability proportional to M 2 Q 2 in the high energy limit (contribu- 
tions of angular velocity and electric potential have been neglected since the black 
hole discharges and finishes its rotation much faster than it evaporates). PBHs 
have been investigated for many different purposes, including tests for quantum 
gravity [2] that are especially active nowadays. 

As was shown by MacGibbon and Webber [3], when the black hole tempera- 
ture is greater than the quantum chromodynamics confinement scale Aqcd, quark 
and gluon jets are emitted instead of composite hadrons. This should be taken 
into account when computing the cosmic-ray flux expected from their evapora- 
tion. Among all the emitted particles, two species are especially interesting : 
gamma-rays around 100 MeV because the Universe is very transparent to those 
wavelengths and because the flux from PBHs becomes softer (oc E~ 3 instead of 
oc E~ v ) above this energy, and antiprotons around 0.1-1 GeV [4] because the 
natural background due to spallation of protons and helium nuclei on the inter- 
stellar medium is very small and fairly well known. This article aims at taking 
into account the contribution from blazars as well as from normal galaxies in the 
gamma-ray background to reduce the available window for PBHs. 

Computing the contributions both from the direct electromagnetic emission 
and from the major component resulting from the decay of neutral pions, the 
gamma-ray spectrum from a given distribution of PBHs can be compared with 
measurements. 

The flux on Earth can be written as: 



dEdt 2 Jt form \R form J JM,{t)dEdV R dMidV 

where tf orm is the formation time, t is the age of the Universe, r is the optical 
depth, R(t) is the scale factor of the Universe at time t, (f> is the individual gamma 
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spectrum from a PBH and d 2 n/dMidV is the initial mass spectrum. When com- 
pared with the observations, this translates (in units of critical density) into [5] : 
Qpbh(M*) < 1.0 x 1CT 8 , bettering substantially previous estimates [6]. 

This limit can be improved when taking into account the "guaranteed" gamma- 
ray background. Computing the contribution from unresolved blazars and the 
emission from normal galaxies, Pavlidou & Fields [7] have estimated the min- 
imum amount of extragalactic gamma-rays which should be expected. The 
first one was computed using the Stecker-Salamon model and the second one 
is assumed to be proportional to the massive star formation rate (which is it- 
self proportional to the supernovae explosion rate) as it is due to cosmic-ray 
interactions with diffuse gas. This background at 100 MeV is estimated at 
&th — 5.45 x 10~ 14 cm _3 GeV _1 . Using the Carr & MacGibbon estimation at 
the same energy [5] §pbh = 7.5 x 10~ 6 Q PBH cm _3 GeV _1 , a new limit can be ob- 
tained by requiring that §pbh + &th < & egret where $ egret is the measured 
flux [8] . To evaluate this later in a very conservative way, both the normalisation 
and the spectral index were chosen (within the error bars) in this paper at the 
value leading to the highest 100 MeV flux : $ egret < 7.94 x 10~ 14 cm _3 GeV _1 . 
This leads to n PBH (M*) < 3.3 x 10~ 9 . 

From the cosmological point of view, this new limit improves directly the 
estimates of the maximum allowed PBH mass fraction (3: 

I C&max (Tu(tl) V'" 

P(M H ) = -=± / e ^> d5 » I} k) e 5 W , 
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where is the horizon crossing time for the considered mode, 5 the density 
contrast (5 m i n ~ 0.7), M# is the Hubble mass at tk and cr'jjitk) = cr 2 (R)\t k 
where <J 2 (R) =< (^f)jj > is computed with a filtering window function with 
R = ^—\t k - The latest computations [9] relate this value to the density parameter 
by 



n PBH (M)h 2 = 6.35 x 10 16 x p(M) 



M 



Our new limit leads to /3(M*) < 1.3 x 10 -26 which is compatible with antiprotons 
estimations [10] and remains the only observational access to such small scales in 
the early Universe. 

Measurements from the GLAST satellite and more refined theoretical predic- 
tions on the background could slightly improve these results but both gamma- 
rays and antiprotons seem to have closed their detection windows. One of the last 
hopes could reside in antideuterons which are very rarely induced by spallation 
below 1 GeV for kinematical reasons [11]. 
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